Electron spin resonance in S = 1/2 antiferromagnets at high temperature 
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We study the electron spin resonance (ESR) of low-dimensional spin systems at high temperature, 
and test the Kubo-Tomita theory of exchange narrowing. In finite-size systems (molecular magnets), 
we found a double-peak resonance which strongly differs from the usual Lorentzian. For infinite 
systems, we have predictions for the linewidth and lineshape as a function of the anisotropy strength. 
For this, we have used an interpolation between a non-perturbative calculation of the memory 
function at short times (exact diagonalization) and the hydrodynamic spin-diffusion at long times. 
We show that the Dzyaloshinskii-Moriya anisotropies generally induce a much larger linewidth than 
the exchange anisotropies in two dimensions, contrary to the one-dimensional case. 



INTRODUCTION 



The paramagnetic resonance is a well-known phe- 
nomenon resulting from the collective precession of the 
total magnetization about an external magnetic field. 
For spin systems with interactions of the standard SU(2) 
Heisenberg form, the total magnetization is conserved 
and there is no relaxation of the magnetization, regard- 
less of the strength of the interaction. It is no longer 
true when anisotropic interactions are present, and the 
magnetization relaxes with a characteristic time scale. 
A theoretical issue is to relate this time scale (or the 
width of the resonance) to the anisotropy strength in 
a many-body system. In strongly interacting systems, 
small anisotropies (as in transition metal compounds) 
lead to an exchange-narrowed Lorentzian resonance^ 
The linewidth is given, generically, by the theory of 
Kubo-Tomita^ and cast into a more general formalism by 
Mori 3 - and Zwanzigi^ However, the Markovian assump- 
tion used by Kubo-Tomita was later argued to break 
down in low dimensional systems. For one dimensional 
systems, deviations from the Lorentzian lineshape were 
indeed observed experimentally and attributed to spin 
diffusion;^ which was not taken into account in the first 
place. Recently, Oshikawa and Affleck^ have developed 
a direct approach in one-dimensional S — 1/2 systems, 
based on an effective field-theory, which is valid at low 
temperatures and could account successfully for the (low- 
)tcmpcrature dependence of the measured lincwidths. It 
is, however, clear that a theory of the high-temperature 
linewidth would be very useful to extract the anisotropies 
experimentally. The Kubo-Tomita formula has been 
widely used in this regime but not thoroughly tested. 

It is nonetheless known that, even in the high- 
temperature regime, a direct application of the Kubo- 
Tomita formula fails in the case of Dzyaloshinskii-Moriya 
anisotropy in one dimension^^ which is first-order in 
the spin-orbit coupling in S — 1/2 systems. It predicts 
indeed a linewidth varying like D 2 /J, which was orig- 
inally argued to be responsible for the large linewidth 
observed in particular in CuCeC ^ 10 ' 11 It is in fact in- 



compatible with a general argument in one dimension, 
which leads to predict a smaller D 4 /J 3 ££ It is then 
of the same order of magnitude as the contribution 
of the exchange anisotropy (as also happens for other 
observables^) , i.e. fourth-order in the spin-orbit cou- 
pling and, hence, small. As alternatives, unconvention- 
nal superexchangei 3 - or dynamical Dzyaloshinskii-Moriya 
interactions^ were invoked to explain the strong broad- 
ening. Nonetheless, it is expected, on general grounds, 
that one-dimensional systems may possess an inherent 
strong broadening because of spin-diffusion at high tem- 
peratures. One issue is to obtain a reliable quantita- 
tive estimation of the effect. Approximations, such as 
the random phase approximation (RPA)^ perturbative 
expansions^ or self-consistent RPA^ suggested in par- 
ticular departures from the Lorentzian lineshape, but it 
is difficult to judge how quantitative these theories are. 

In contrast, in two spatial dimensions, it is surpris- 
ing to notice (empirically) that applying the Kubo- 
Tomita's formula seems to give rather accurate results 
at high temperatures. For instance, in the Shastry- 
Suthcrland compound, SrCu2 (B0a)2, this gives an out- 
of-planc Dzyaloshinskii-Moriya interaction D = 2.4 K,— 
while neutron scattering gives 2.1 KJ^ 

In the following, we will give a quantitative estima- 
tion of the linewidth in one and two dimensions and the 
corresponding lineshapes. We predict an enhancement 
of the linewidth with respect to the perturbative Kubo- 
Tomita formula, particularly strong when the diffusive 
motion is assumed at long times. Moreover, while we 
agree that the linewidth induced by the Dzyaloshinskii- 
Moriya anisotropy or the exchange anisotropy are of the 
same order in one dimension (or more generally for the 
reducible^ Dzyaloshinskii-Moriya interactions) , we claim 
that, in two dimensions, the irreducible Dzyaloshinskii- 
Moriya interactions lead to a linewidth essentially vary- 
ing like D 2 / J, i.e. like in Kubo-Tomita, with a prefactor 
that we shall estimate. 

An approach consists of exact numerical diagonaliza- 
tions of the Hamiltonian and a calculation of the dynam- 
ics using the Kubo formulai 2 ^ The approach is interesting 
not only because there is no approximation made but also 
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because it provides information on the whole ESR spec- 
trum and tells us whether there are more resonances than 
the paramagnetic resonance. It is of course limited by the 
system size and especially at finite temperatures where 
one needs to calculate all eigenstates. As we shall see, 
a direct computation of the Kubo formula gives results 
which arc difficult to extrapolate to the thermodynamic 
limit. Since, however, the calculation is exact for small 
sizes, the prediction for the lineshape may be interesting 
for molecular magnets, in particular to disentangle dif- 
ferent mechanisms. However, to obtain information for 
bulk systems, we suggest that it is more interesting to 
calculate the memory function by exact diagonalization, 
in the Mori-Zwanzig framework. While the calculation 
of the memory function is exact only at short-times for 
finite-size systems, its validity goes beyond the perturba- 
tive short-time expansion developed earlier.— ^ We find 
a clear slowing down of the memory function compared 
to Kubo-Tomita in one dimension (but not in two di- 
mensions) , possibly indicative of a crossover to spin dif- 
fusion. We can then test the spin diffusion assumption 
for the long-time behaviour. It is indeed reasonable to 
assume that the short-time behaviour depends on the mi- 
croscopic details of the lattice (and hence needs a precise 
calculation) while the long-time may be more universal 
(and depends primarily on the spatial dimension). This 
allows to obtain a quantitative idea of how accurate the 
Kubo-Tomita's result is in one and two dimensions. 



II. KUBO FORMULA AND FINITE-SIZE 
SYSTEMS 

The linear response of the spin system to a long wave- 
length oscillating magnetic field along, say, an cc-axis is 
given by the Kubo formula (or susceptibility), 

r+oo 

X (u) = i cbje^([S x (t),S x (Q)}) (1) 
Jo 

where S x = Sf is the total magnetization along the x- 
axis. The absorption cross-section of the incident electro- 
magnetic wave is proportional to wx"(w), and 

x»^(i- e -^)s*», (2) 

where 

/+oo 
du)e iut {S x (t)S x (0)) 
-OO 

= \ E e-^\(n\S x \p)\ 2 S(u; -E p + E n )(3) 

with a sum over the eigenstates \n) and \p) of the Hamil- 
tonian {Z is the partition function). At high tempera- 
tures (T > J), = (uj 2 /T)S xx {oj) and we show 
only S xx (ui) in the following. Furthermore, all the states 
have the same weight (e~^ En ~ 1) and S xx (lu) is an 



even function of ui (S xx (oj) = S xx (—lj)). To com- 
pute Eq. ©, we have calculated all eigenenergies and 
eigenstates of the Hamiltonian by exact numerical di- 
agonalization. We have first considered the Heisenberg 
Hamiltonian for one dimensional spin rings (with periodic 
boundary conditions) with two types of SU(2)-symmetry- 
breaking anisotropics, one is the Dzyaloshinskii-Moriya 
interaction (first-order in spin-orbit coupling), 

H = ^JS l .S l+1 +D. i .(S l x S i+1 ), (4) 

i 

where is a staggered vector from bond to bond 
(strenght D) and taken along the z-axis. The second 
is the XXZ exchange anisotropy (second-order in spin- 
orbit coupling), 

W = J^-^+i - tJS?S* +1) (5) 

i 

with an easy-plane anisotropy (S > 0). Importantly, H 
can be mapped onto W by rotating the spin operators, 
5+ = e^-^Sf with 9 = D/2J, giving 5 = D 2 /2J 2 at 
the lowest order i 8 ' 9 ' 12 In this mapping, the long wave- 
length oscillating magnetic field acquires a staggered 
component. In the following, we restrict ourselves to an 
external field h along z. In this case, the z-component of 
the total magnetization, S z , is conserved and the mag- 
netic field splits the excited states according to —hS z but 
does not change the matrix elements in Eq. ([3]) . In addi- 
tion, the operator S x changes S z by ±1, so that E p —E n is 
shifted by ±h. As soon as h is larger than the linewidth, 
the h 7^ lineshape is the same as for h = shifted by 
±h. 

The result of the calculation of the Kubo formula §3§ 
for N = 16 spins interacting with the Hamiltonian (U) 
is shown in Fig. Q] (the delta peaks in ([3]) arc slightly 
broadened by a Gaussian formal). We find a broad res- 
onance centered at the magnetic field h (there is no fi- 
nite shift). Since we are at high temperatures, we stress 
that S xx (u>) is the same around u = —h. The calcu- 
lation being exact, this is the lineshape we predict for 
molecular magnets of the ring geometry, at high tempera- 
ture, would the mechanism giving the linewidth be purely 
magnetic in origin. In Fig. [U the different curves cor- 
respond to different strengths of Dzyaloshinskii-Moriya 
coupling. They are normalized in intensity and rescaled 
in frequency using the linewidth, noted cjq, and defined 
as the half-width at half-maximum of the full line. They 
all collapse onto a single (quasi)-universal lineshape with 
a double-peak structure, in a wide range of anisotropy 
strengths, D ~ 0.1 — 0.8 J and sizes. If we consider the 
special case h = 0, the double-peak structure can be seen 
as a pseudo-gap. The prediction of this lineshape seems 
at odds with the conventional exchange-narrowing the- 
ory which, in general, predicts a single Lorentzian. We 
will discuss the origin of these differences below. In addi- 
tion to the double-peak structure, there is an additional 
weight in the wing. In fact, the line-shape can be fitted 
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FIG. 1: (Color online) ESR lineshape of a one-dimensional 
chain of N = 16 spins (T J). The curves are rescaled with 
the half- width at half-maximum, u>o (inset). 

by four Lorcntzians, two centered at ±0.59o;o and two at 
±3(0.59o;o) (each component is given by a dashed line in 
Fig. [T]and the fit by a solid thick line). 

The linewidth ojq increases basically like aD 2 / J (inset 
of Fig.[T]) with a rather small prefactor, a. The prcfactor 
is size-dependent and is difficult to extrapolate to the 
thermodynamic limit. For this reason (and following the 
discussion of section llllj) , the special double-peak line- 
shape may or may not survive in the thermodynamic 
limit. For molecular rings, this should be relevant and 
the prefactor is found to be a = 0.090 (N = 12), 0.071 
(N = 14) and 0.059 (N = 16). 

We now emphasize the difference between the 
Dzyaloshinskii-Moriya model and the XXZ model, as far 
as the response is concerned. In Fig. [3J we show the res- 
onance line in absolute values for the two models with 
coupling strengths related by the mapping 28 = (D/J) 2 . 
The two resonance lines are undistinguishable at low en- 
ergies (inset of Fig. [2]) but differ at high energy in the 
background (Fig. [2]). The difference comes from the 
small oscillating staggered field (or umklapp contribution 
from the spin fluctuations at q = ir) that the mapping 
induces^ 

S xx («) = cos 2 9S XX («) + sin 2 6S XX (w) (6) 

where S xx and S xx correspond to the Dzyaloshinskii- 
Moriya and XXZ model, respectively, and S xx (uj) is the 
q = 7r response. The second term gives extra "forbidden" 
resonances (or "satellite" lines) at low temperatures.— 
These resonances have disappeared at high temperatures. 
Instead, this contribution gives a broad background, cov- 
ering the whole energy range. This is the origin of the dif- 
ference of the two curves at high energy. To summarize, 
the resonance lineshape as shown in Fig. [T] is the same 
for both the Dzyaloshinskii-Moriya and the XXZ models 
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FIG. 2: (Color online) ESR lineshape (high-energy tail) for a 
one-dimensional chain (T ^> J) for the Dzyaloshinskii-Moriya 
and XXZ models (using 28 = D 2 /J 2 ). The response is iden- 
tical only at low-energy (inset). 

as long as we use the correspondence 25 = (D/J) 2 (and 
not too small intensities). 

Note that we have found deviations about the univer- 
sality of the low-energy lineshape at strong anisotropics 
(cither Dzyaloshinskii-Moriya or XXZ). This is expected 
because, for 5 = 1, the model is XY and the lineshape 
is known to be exactly Gaussian^ 2 - We have recovered 
this result numerically: we have found no difference with 
a Gaussian at 5 = 1 (no finite-size effects). The way it 
crossovers from a Gaussian function to four Lorentzians 
(by reducing 5 from 1) is to depress the response at lo = h 
and let side peaks appear, which decrease in height. It is 
therefore clear that in the range of strong anisotropics, 
the lineshape is not universal at all. 

III. KUBO-TOMITA THEORY, EXACT 
CALCULATION OF THE MEMORY FUNCTION 

We now discuss the Kubo-Tomita theory, in order to 
understand why the lineshape is not a single Lorentzian 
in such finite-size systems and go beyond the limita- 
tions imposed by finite-size calculations. A way to clarify 
the assumptions is to start with the more general Mori- 
Zwanzig formalism*^ In this framework, the spin fluc- 
tuations (those that do not depend on time through the 
magnetization function) are viewed as a random noise 
(the bath) that forces the magnetization to equilibrate. 
The system is then governed by a Langevin equation, 

j t S x {t) = -J* dr£(i - r)S x (r) + rj(t) = F(t), (7) 

where the "memory function" (or self-energy) £(f) is re- 
lated to the correlation function of the random force rj(t) 
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by the fluctuation-dissipation theorem. However, neither 
E(i) nor rj(t) are easily expressed in terms of the orig- 
inal spin Hamiltonian because of unknown projection 
operators It is only at second-order in perturbation 
theory (here in the anisotropy) j» that E(i) is given by 
the total force correlation function 

E(t) = (F(t)Ft(O)) (8) 

where the total force (or torque) is given by F(t) — 
ie %ut [H, S x \e~ l ' Ht . Since F(t) is linear in the anisotropy, 
E(i) is second-order, and perturbation theory holds for 
small anisotropics. Eq. © is precisely what Kubo and 
Tomita have obtained directly by perturbation theory. 

The second assumption is that S x (t) is slow compared 
to the relaxation time of the force (local equilibrium)^ 
It is true that the force F is not conserved by the dy- 
namics (whereas S x (t) is nearly conserved) and evolves 
on a short time-scale of order \ j J but the absence of 
memory at long times is an assumption. If it is true (the 
system is Markovian) , Eq. (J7J simply defines a relaxation 
time, t c (or Onsager's transport coefficient) 



1 f +ao 

- = / drE(r). 

r c JO 



(9) 



In this "fast" regime, the lineshape is therefore always 
Lorentzian with half-width at half-maximum (lincwidth) 
given by l/r c . In contrast, if slow hydrodynamic modes 
are present, the spin motion may be diffusive^'^ In this 
case, the memory function has a power-law tail E(t) ~ 
\jt d l 2 and Eq. fl5J diverges for d < 2. In fact, this signals 
a change of lineshape and the general solution of for 
the correlation function is 



S xx (uj) 



1 



E» 



4tt [u + £'H] 2 + E"(w) : 



(10) 



with E(w) = i J Q E(t)e IW *<ii the Laplace transform of 
the memory function, and E'(w), E"(cj) are its real and 
imaginary parts. If "slow" processes are present, E(cj) 
varies considerably near u> = and the lineshape devi- 
ates from a Lorentzian. The form of E(£) is therefore 
of importance and various assumptions have been used. 
The original (third) assumption of Kubo- Tomita was to 
use a Gaussian decay for E(f), based on the short-time 
evolution given by perturbation theory 



m = E 



(itr 



M, 



n+2 



(11) 



where M n+ 2 are the moments, M 2 = E(0) = (F 2 ), 
M 3 = ([%, F]F), M 4 = ([U,F][F,H]) etc.. They 
can be calculated exactly at high temperatures. For 
the XXZ model ©, M 2 /J 2 = <5 2 /4, M 3 = 0, and 
M4/J 4 = 3<5 2 /8 - 8 3 /4 + 5 4 /4:. So at short times, 
E(i)/E(0) = 1 - iw e 2 t 2 with the definition w 2 = 
Kubo and Tomita have postulated that 
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FIG. 3: (Color online) Memory function £(t) at short times 
(N is the system size, 5 = 0.08). Convergence to N — 00 is 
obtained for tJ < 6.5. The Kubo- Tomita gaussian holds for 
tJ < 1.5. The interpolation to the long-time spin-diffusion 
1/t 1 / 2 is shown (dashed line). 



with ui 2 matching that of the short-time expansion (|11[) 
at second order, uj 2 = M 4 /M 2 = 3 J 2 / 2 - 5 J 2 + 
S 2 J 2 . Eq. (fT2"j) then basically corresponds to assuming 
a ressummation of an infinite number of terms. The 
Kubo- Tomita formula for the linewidth is then at high 
temperature, 



1 

Tc. 



7TM 2 

2 oj k 



nS 2 
3 4 



-J 



(13) 



To summarize, the assumptions that lead to Eq. (|13[) 
are (i) perturbation theory in 8 <Si 1, (ii) Markovian be- 
haviour, and (iii) perturbation theory in time ( C 1/J. 

In order to go beyond perturbation theory and test 
these assumptions, we have calculated 



s(*) = ^E e "^IH[« ) 5 B ]|p>l a c < ( B " 



-E p )t 



(14) 



II . p 



E(i) = E(0)exp(-a; 2 i 2 /2) 



(12) 



by exact diagonalization of the Hamiltonian [sec Eq. ((3|) 
for the definitions]. We restrict ourselves again to high 
temperatures (e~^ En = 1) and have computed all matrix 
elements. Here we calculate quantities for the model (|5j) 
which makes no difference as we mentioned. At short 
times (Fig. [3]), E(i) decreases indeed on a time scale 
~ 1/J and the Kubo- Tomita gaussian decay holds for 
tJ < 1.5 (pcrturbativc regime). For tJ > 1.5, E(f) de- 
creases slowly and by comparing the results for differ- 
ent sizes, we see no difference between the N = 16 and 
N = 18 curves for tJ < 6.5. We can therefore be con- 
fident that for tJ < 6.5, E(i) is representative of the 
thermodynamic limit. The result is different from the 
gaussian decay and there is a considerable slowing down 
for 1.5 < t < 6.5. We have here access to a regime that 
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would be difficult to access by perturbation theory. We 
will come back on this non-perturbative regime in the 
next section. Above tJ > 6.5, finite-size effects domi- 
nate. In particular at longer times (Fig. 2]), £(£) acquires 
an additional broad peak at about t a J ~ 10 — 17. t a 
increases almost linearly with the length of the chain N, 
and its amplitude decreases with ~ 1/iV 3 / 2 (as shown 
by the fit in the inset of Fig. 0]). A possible inter- 
pretation is that this additional peak is related to the 
torque travelling around the ring. In this case, indeed 
the arrival time is expected to be controlled by J and 
weakly sensitive to the anisotropy, which is what we 
found. In fact if we take the des Cloizeaux-Pearson 
speed of the hydrodynamic excitation, the arrival time 
is t a = N/c = 2N/ir ~ 6 — 11. This additional peak to- 
gether with the long-time decay should be at the origin of 
the "double-peak" structure observed in Fig.[TJ A simple 
way to see the effect of a long-time decay of S(i) is to 
consider the gaussian memory function of Eq. (|12p with 
a long-time scale l/uj e ^> 1/J. The Laplace transform, 
is a sum of a Gaussian function (imaginary part) 
and a Dawson function (real part), and S xx (u>) [Eq. (ITU1) ] 
acquires a double peak structured This occurs when the 
time scale of the memory function becomes comparable 
to that of the spin motion l/uj e > t c \/2tt/(2 + y/n). For 
a general £(t), the change of sign of the curvature of 
Eq. (fTU| at u) = signals the occurence of a double- 
peak structure. This happens when r c < (t) + (i 2 ) 1 / 2 
where (t n ) = (1/n!) / +o ° E(t)t n dt/ f+°° £(t)tft. As a 
consequence, there may be a critical size above which the 
double-peak structure disappears. These long-memory 
processes (non-Markovian) explain why deviations from 
a Lorcntzian lineshapc take place in finite-size systems. 
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FIG. 4: (Color online) Memory function E(t) at long times 
(system size N, S = 0.08). There is a peak at t a J ~ 10 — 17 
whose amplitude decreases with N (inset). The assumed long- 
time spin-diffusion tail is also shown (dashed line). 



IV. RESONANCE IN ONE DIMENSIONAL 
THERMODYNAMIC SYSTEMS 

In the previous sections, we have given the exact line- 
shape of a finite-size system. As we have explained, it 
is difficult to extrapolate it in the thermodynamic limit. 
We now discuss a way to obtain the resonance line in the 
thermodynamic limit by using some additional assump- 
tions. The crucial point is to obtain the memory function 
£(i), but it cannot be computed at any time t by exact 
diagonalization. It is therefore interesting to assume a 
long-time behavior and test its consequences on the reso- 
nance line, which can then be compared to experiments. 
For this, we will interpolate between the exact result at 
short times (up to the time where £(i) has converged 
for N — > oo, i.e. t J — 6.5 in Fig. [3j and the assumed 
hydrodynamic behavior at long times. 

First we note that the slowing down obtained in the 
non-perturbative regime (1.5 < tJ < 6.5) already leads 
to increase the linewidth compared with Kubo-Tomita. 
If we assume for instance an abrupt cutoff at t J = 6.5 (re- 
member that for t J > 6.5, S(t) is dominated by finite-size 
effects), Eq. ([9]) defines a lower bound for the linewidth. 
The lower bound, fitted by ~ 0.41<S 2 , is yet larger than 
the Kubo-Tomita formula, 0.26J 2 [Eq. ([15])] (both shown 
in Fig. P. 

We now consider the idea that the long-time be- 
havior may be universal and governed by a spin dif- 
fusion equation. We recall that the diffusion of £(t) 
is obtained assuming a RPA decoupling of the 4- 
spin correlations in Y,{t)p& {S+{t)S*{t)S^{Q)S^{fi)) ~ 
(S+{t)S~,{Q))(S z q (t)S z ,(0)). Since S z is conserved, 
the long wavelength modes are supposed to be dif- 
fusive with (S z {t)St q (0)) ~ e^A On the other 
hand, (S+(t)Sl q (0)) ~ e -o^t e -t/T B is cut _ off by the 
anisotropy, where t c is precisely the characteristic time- 
scale we are looking for, as emphasized by Rciter and 
Boucher^ When summed over q, this leads at long times 
(t > 1/J) to 

p-t/ro 

S W - (15) 

where To is the cut-off time of the memory function: 
within the self-consistent RPA, tq = t c . As explained 
above, we cannot test this behavior quantitatively by 
exact diagonalization, although the slowing down for 
1.5 < t J < 6.5 seems to indicate a crossover regime. To 
test the idea of spin-diffusion, Fabricius and McCoy have 
examined directly the two-spin autocorrelation function 
for the N = 16 chain at high temperatures and con- 
cluded that the exponent may be closer to 0.7 in one 
dimension and increases with ^ So the exponent 1/2 
may be underestimated, thus leading to overestimate the 
linewidth. Since r c now explicitely enters in Eq. (fl~5"j) . 
Eq. © is a self-consistent equation. This is precisely 
the equation that Reiter and Boucher solved^ Instead 
of relying on RPA to obtain the prefactors, here we shall 
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FIG. 5: (Color online) Linewidth vs. anisotropy for infinite 
one-dimensional chains (lower bound, Kubo-Tomita and self- 
consistent calculation) . 



FIG. 6: (Color online) ESR lineshape for infinite one- 
dimensional chains, assuming different cutoffs to of the dif- 
fusive tail. The arrow (inset) indicates the self-consistent t c . 



assume a long time tail of the form (TT5l) and interpolate 
to the short-time behavior we have obtained by exact di- 
agonalization. The idea of interpolating the two types of 
behavior was used by Gulley et ai, together with per- 
turbativc short-time expansion.— The interpolation pro- 
vides a "variational" function E(i,T c ) and we solve (in 
a similar spirit as in Ref. fltil ) the (now) self-consistent 
Eq. (|9]) with the new E(i,r c ). In fact the long-tailed 
£(i) changes the lineshape according to Eq. (TlO|) . In this 
case, Eq. (JU) does not hold anymore and we calculate 
the linewidth using instead the self-consistent equation, 
S xx (1/t c ) = S xx (0)/2 where S xx (uj) is the solution given 
by Eq. (|10|) . The difference between the two procedures 
is in fact quite small (from 4% to 7% depending on 5), 
and we show only the full self-consistent result in Fig. [5] 
(squares). This provides an estimation of the linewidth 
for a pure one dimensional model within the spin diffu- 
sion assumption. 

It is now possible to relax the assumption that the 
cutoff To = r c in Eq. (Tj"5|) . It does not seem to be phys- 
ically reasonable to take To > r c at high temperatures 
because it would mean having finite torque-torque corre- 
lations but no spin-spin correlation after a certain timed 
In this case, Eq. (|9]) leads to a divergence of the linewidth 
~ y^To whereas the use of the correct lineshape provides 
a less diverging result (shown in the inset of Fig. [5]). The 
other assumption tq < t c is physically more relevant for 
two reasons. First tq may be an intrinsic characteris- 
tic time of the one-dimensional model that terminates 
the spin-diffusion before the anisotropy cut-off: To = t c 
is valid only at the RPA level. Second, r may be an 
extrinsic "noise" : there are additional interactions, such 
as interchain couplings, that also tend to terminate the 
diffusioni^r— In fact it was argued that tqJ depends 
on the interchain coupling J' through (J 7 / J)~ 4 / 3 , which 



may indeed be a shorter time-scaled On the other hand, 
following the exact result at short times, the cutoff can- 
not be chosen smaller than tqJ = 6.5 (but of course if 
J'/ J becomes very strong, then one has to recalculate 
the short-time behavior as well) . The linewidth decreases 
when to decreases; as a consequence, we regard the re- 
gion represented in Fig. [5] by a hatched area as possibly 
relevant, experimentally. As we now explain, in this re- 
gion, the lineshape is a function of To, thus giving the 
possibility to say how long the motion is diffusive and if 
it is diffusive at all. 

In Fig. [6l we give the line-shapes (normalized by the 
area), computed from Eq. (Ti"0|) . for different cut-offs of 
spin-diffusion, To. At small To <C t c , the lineshape is 
always Lorcntzian (straight dashed line in the inverse 
representation on the right). For the self-consistent 
cut-off to = t c , the lineshape strongly deviates from a 
Lorcntzian. The deviation is less than that found by us- 
ing the RPA calculation with a spin-diffusion tail [which 
is given by the Fourier transform of exp(— At 3 / 2 ) (see in- 
termediate dashed line)]£ and still far from the gaussian 
profile (upper dashed line). For t > t c we find a "double 
peak" structure similar to that found in finite-size chains. 
We have therefore two possibilities for the double-peak 
structure as a function of the system-size: either it sur- 
vives in the thermodynamic limit (which means that for 
some reasons To could remain larger than t c ) or it dis- 
appears. For a general To, the inverse lineshape changes 
smoothly as function of t and therefore a comparison of 
the lineshape with experiment would provide a measure 
of the cutoff time, to- 

The linewidth (Fig. [5]), together with the lineshape 
(Fig. [5]) , can therefore be used to extract the anisotropy 
strength and test the spin diffusion assumption. Note 
that in one dimension, because of the mapping we 
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FIG. 7: (Color online) Memory function for the two- 
dimensional S = 1/2 kagome antiferromagnet: effect of the 
Dzyaloshinskii-Moriya strength (left) and finite-size effects 
(right). 



have discussed, it is not possible to tell whether the 
lincwidth is due to the Dzyaloshinskii-Moriya or ex- 
change anisotropics both contribute to an equal 
amount, using the replacement 28 = (D/J) 2 . 



right). If we had only the N = 12 cluster, it would 
be tempting to conclude that Kubo-Tomita is essentially 
exact in two dimensions. However, the N = 15 cluster 
deviates from Kubo-Tomita. It is not clear whether this 
is an artifact due to the odd size of the cluster. We now 
proceed with the interpolation to the spin diffusion as- 
sumption (which is in 1/t in two dimensions). Since we 
see no clear slowing down in the present case, we have 
used two different times for the interpolation, tJ = 1 
(which certainly gives an upper bound) or t J = 1.5 (see 
the two dashed lines in Fig. [7]). 



0—0 Self-consistent calculation, 1st inlerpol. 

HI— M Self-consistent calculation, 2nd inlerpol. 
- - - 1/x =0.32 (D/J) 2 | 0.22(D/J) 2 

Kubo-Tomita: 1/T =0.1 1(D/J) 2 

_ - ZnCu,(OH) 6 Cl 3 , Zorko el al. 




V. RESONANCE IN TWO DIMENSIONAL 
THERMODYNAMIC SYSTEMS 

In two dimensions, the pattern of Dzyaloshinskii- 
Moriya vectors may generally be more complicated, de- 
pending on the symmetries of the crystal structure. In 
particular it is not always possible to perform the map- 
ping onto an exchange anisotropy because of the frus- 
tration that the closed loops of the lattice introduce (ir- 
reducible components). In this case, there seems to be 
no reason why the Dzyaloshinskii-Moriya term should 
contribute to the linewidth by an equal amount as the 
exchange anisotropy, and we shall see that it does not. 

We have considered the example of the kagome anti- 
ferromagnet where the pattern of Dzyaloshinskii-Moriya 
vectors precisely have irreducible components^ These 
are the z-component of the Dzyaloshinskii-Moriya field. 
We have calculated the memory function E(f) by exact 
diagonalizations of clusters of up to 15 spins. We see 
that the result is weakly dependent on the strength of 
the Dzyaloshinskii-Moriya coupling, D (Fig. [71 left), ex- 
cept through S(0) = D 2 /8 (at high temperatures). As 
a consequence, the linewidth should vary essentially like 
D 2 /J. We now proceed like in one dimension to obtain 
a quantitative estimation of the prefactor. In fact finite- 
size effects are somehow more difficult to handle because 
we have less sizes available. In addition, we see a par- 
ity effect between clusters of odd and even sizes (Fig. [7J 



FIG. 8: (Color online) Linewidth vs. anisotropy for the 
S = 1/2 kagome antiferromagnet (Kubo-Tomita and self- 
consistent calculation). The horizontal line shows the exper- 
imental linewidth of ZnCu 3 (OH) 6 Cl 3 (Ref.[H). 
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FIG. 9: (Color online) ESR lineshape for the two-dimensional 
kagome antiferromagnet, assuming different cutoffs to of the 
1/t diffusive tail 
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The two results of the self-consistent calculation are 
shown in Fig. [5] (squares) and are fitted by l/r c = 
0.32( J D/J) 2 and 0.22(D/J) 2 for small anisotropy. The 
linewidth is therefore larger than in one dimension. To 
put numbers, if we take D = 0.2J, the linewidth is 
6.0 x 10 -3 J (Fig. ©, whereas for the same D in one di- 
mension, we have to take S = (D/J) 2 /2 = 0.02, and the 
linewidth is 1.0 x 10~ 3 J (Fig.©. Now, as in one dimen- 
sion, the area below the self-consistent result in Fig. [8] 
could be experimentally relevant (because of interplane 
couplings, for instance). Away from the Kubo-Tomita 
line (lower bound in Fig.©, we can see in Fig.^that the 
lineshape starts to deviate from a Lorentzian (but less 
than in one dimension). 

We now apply this to the 5 = 1/2 kagome compound 
ZnCu3(OH)6Cl3.— The analysis of the linewidth using 
the Kubo-Tomita formula provided D/J = 0.08^ If we 
assume that spin-diffusion takes place, the self-consistent 
calculation reproduces the same experimental linewidth 
providing 0.044 < D/J < 0.08 (see dashed line in Fig.©). 
It is difficult to compare directly with the ESR lineshape 
of ZnCii3(OH)gCl3 because the experiment was done on 
powder samples.— One needs to average over all direc- 
tions of the field and we have assumed a single direc- 
tion here. Nonetheless, Fig. [TU] shows that the devia- 
tions from a Lorentzian in the wings (here we restrict 
to w < ft, which is not spoiled by an impurity line) can 
be accounted for by the spin diffusion assumption and 
D/J^ 0.05. The result may be partly coincidental be- 
cause there are other sources of deviations, such as the 
chemical disorder and the polycristalline nature of the 
sample. Still we note that this coupling strength is not 
incompatible with the analysis of NMR which stated that 
D j J ' > 0.05.— In any case, from the ESR linewidth alone, 
we can conclude that 0.044 < D/J < 0.08. 



i — r 




CO-h (Tesla) 



FIG. 10: (Color online) ESR signal (derivative of the reso- 
nance line at to < h): theory and experiment on the kagome 
S = 1/2 ZnCu 3 (OH 3 )Cl 2 (from Ref.tH). 



VI. CONCLUSION 

First, for finite-size chains (0 dimension), we have 
found a special (and exact) "double-peak" lineshape 
(Fig. © which may be interesting in molecular magnets 
having the geometry of a ring, especially because the line- 
shape seems insensitive to the strength anisotropy (up to 
a rescaling) and size (in the limit of small size). We have 
stressed that such a lineshape usually results from non- 
Markovian correlations in the memory function, which 
we did find in the calculation of £(i) in this case. 

In the thermodynamic limit, we have obtained the elec- 
tron spin resonance in one and two dimensions for the 
anisotropics relevant to S = 1/2 systems, by comput- 
ing the memory function. We have provided the exact 
numerical result for £(t) at short times and used it as 
an initial condition for a long-time spin-diffusion tail, 
which was taken as an assumption. We have then cal- 
culated the characteristic relaxation time (or linewidth) 
self-consistently, as well as the lineshape. 

In one dimension, the spin-diffusion assumption leads 
to enhance the linewidth by a factor that we have quan- 
titatively estimated (see Fig. [5] [squares]). It is about 
4 times larger than the perturbative Kubo-Tomita's re- 
sult. When this is true, the lineshape strongly differs 
from a Lorentzian. To compare with experiments, it is 
interesting to consider another time-scale, To which is 
the cutoff of the memory function (the physical reason 
may be intrinsic or extrinsic; for instance, because of 
interchain couplings). When To is reduced from its self- 
consistent value, the linewidth decreases and the line- 
shape crossovers to a Lorentzian. A comparison of both 
the linewidth and lineshape should therefore allow ex- 
perimentally to determine the anisotropy strength and 
To (with a better accuracy than the use of the Kubo- 
Tomita formula). Note that, given a resonance line, we 
can extract the anisotropy parameter 5 but we cannot 
tell whether S is due to the exchange anisotropy or the 
Dzyaloshinskii-Moriya field (unless the later is forbidden 
by symmetry): both contribute at the same order of mag- 
nitude in one dimension*^ 

In two dimensions, the situation is different and there 
are two distinct cases: 

• Irreducible case (the Dzyaloshinskii-Moriya vectors 
do not sum to zero while going around closed loops 
of the lattice): we have argued that the linewidth 
should vary essentially like D 2 /J. 

• Reducible case (the Dzyaloshinskii-Moriya vectors 
do sum to zero): an exact transformation maps 
again the model onto an exchange anisotropy and, 
as in one dimension, both contribute at the same 
order (S 2 J~D 4 /J 3 ). 

The result in the irreducible case is in fact what the 
Kubo-Tomita formula gives, qualitatively Quantita- 
tively, the prefactor is bounded below by the Kubo- 
Tomita result. In fact if we restrict to the cluster of 
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even size (N = 12), we see an excellent agreement with 
the gaussian decay and we would be tempted to conclude 
that we see no quantitative difference with Kubo-Tomita. 
By computing the next cluster (TV = 15), however, the 
difference turns out to be larger. While it is difficult to 
know if the difference comes from the odd size of the 
cluster, we have given in the later case the interpolation 
to the \/t d / 2 (d = 2). The self-consistently calculated 
linewidth is in this case close to 0.22D 2 / J or 0.32D 2 / J 
(Fig. [H]) depending on the interpolation time, but still 
larger than the Kubo-Tomita's formula, 0.11D 2 /J (and 
in between, if interplane couplings are present). 

We have applied the present theory to the two- 
dimensional kagome ZnCu3(OH3)Cl2,— and concluded 
that 0.044 < D/J < 0.08 from the ESR linewidth. In 
addition, the spin diffusion assumption seems to account 
for the deviations of the lineshape from a Lorentzian 
(Fig. [TU|) . but single crystals are needed to avoid ad- 
ditional averaging effects. In this compound, the out- 
of-plane Dzyaloshinskii-Moriya component is irreducible 
while the in-plane is reducible^ For single crystals, we 
therefore predict a much larger linewidth when the field 



is perpendicular to the plane. 

We also note that, in absolute values, the irreducible 
Dzyaloshinskii-Moriya interactions in two dimensions in- 
duce a linewidth larger than in one-dimension (where the 
Dzyaloshinskii-Moriya interactions are always reducible) . 
The essential factor that governs the linewidth in low- 
dimensional S — 1/2 systems is not the strength of the 
divergence 1/t 1 / 2 versus 1/t of the spin-diffusion (assum- 
ing it exists) but rather the reducible versus irreducible 
character of the Dzyaloshinskii-Moriya interaction. We 
have illustrated this on the kagome lattice but, at high 
temperatures, the shape of the lattice should not matter 
(as long as there are closed loops with Dzyaloshinskii- 
Moriya vectors not summing to zero). 
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